Curvature-driven coarsening in the two dimensional Potts model 
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We study the geometric properties of polymixtures after a sudden quench in temperature. We 
mimic these systems with the g-states Potts model on a square lattice with and without weak 
quenched disorder, and their evolution with Monte Carlo simulations with non-conserved order 
parameter. We analyze the distribution of hull enclosed areas for different initial conditions and 
compare our results with recent exact and numerical findings for q = 2 (Ising) case. Our results 
demonstrate the memory of the presence or absence of long-range correlations in the initial state 
during the coarsening regime and exhibit super-universality properties. 



I. INTRODUCTION 

Domains formed during the evolution of mixtures are 
of both theoretical and technological importance, appli- 
cations including foams [l|, cellular tissues 0, super- 
conductors 3 1 magnetic domains [1, Q , adsorbed atoms 
on surfaces, etc. In particular, in metallurgy and sur- 
face science, the polycrystalline microstructure, and its 
time evolution, are important in determining the mate- 
rial properties. 

After a sudden change in temperature (or in another 
suitable control parameter) that takes the system from 
the high temperature phase into the coexistence region, 
the system tends to organize in progressively larger or- 
dered structures. The system's temporal evolution is 
ruled by thermal, diffusive, and curvature driven pro- 
cesses, and the actual grovifth \a.w depends on general 
features suchlike the presence of quenched disorder, the 
dimension of the order parameter and whether it is con- 
served or not. Although much is known for binary 
mixtures and systems with twofold ground degeneracy 
{q — 2), much less is understood for polymixtures and 
manifolded ground states {q > 2). In the latter case, 
topological defects pin the domain wall dynamics and 
even in pure systems thermal activation is necessary to 
overcome the corresponding energy barriers. 

Many interesting cellular growth processes are cap- 
tured by a curvature-driven ordering processes in which 
thermal effects play a minor role. These are ruled by 
the Allen-Cahn equation in which the local velocity of 
an interface is proportional to the local curvature, v = 
— (A/27r)/t, where A is a temperature and g-dependent di- 
mensional constant related with the surface tension and 
mobility of a domain wall and n is the local curvature. 
The sign is such that the domain wall curvature is dimin- 
ished along the evolution. In d = 2 the time dependence 
of the area contained within any finite domain interface 



(the hull) on a flat surface is obtained by integrating the 
velocity around the hull and using the Gauss-Bonnet the- 
orem: 
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where Ui are the turning angles of the tangent vector 
to the surface at the n possible vertices or triple junc- 
tions. In some systems, such as the Ising q = 2 model, 
ai — since there are no such vertices and we ob- 
tain dA/dt — —A for all hull-enclosed areas, irrespec- 
tive of their size [1, 0]- In highly anisotropic systems, 
as the Potts model studied here, the angles are all dif- 
ferent. In some of the systems that motivate this work, 
like soap froths, the angles are all equal to 27r/3, that 
is, ai = 7r/3,Vi. Focusing on systems in which there is 
a one-to-one correspondence between vertices and sides 
(meaning that we exclude the Ising limit in which there 
are no vertices but each wall has one side) the above 
equation reduces to the von Neumann law [a, Q for the 
area Ar, of an n-sided hull-enclosed area: 
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n> 1. 
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Whether a cell grows, shrinks or remains with constant 
area depends on its number of sides being, respectively, 
larger than, smaller than or equal to 6. The above law 
can be extended to the case in which the typical internal 
angle depends on the number of sides 0] and to 3d [l3| 
as well. 

Potts models, that is to say q-state spin models on a 
lattice, simulate grain growth. Early studies of the 2d 
q-state Potts model with non-conserved order-parameter 
dynamics have shown that the growth law is analogous 
to that of the (non conserved) Ising model, where in the 
scaling regime the characteristic length scale follows the 
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t^/^ AUen-Cahn law [HI, regardless of the value of q. 
However, for q > 2 and T = 0, this power law growth 
does not hold for long times, the characteristic length 
scale converging to a limiting value. This is in agreement 
with the Lifshitz-Safran criterium ^13i,il4] that states that 
when the ground state degeneracy is q > d + I, where d 
is the system dimensionality, there might be domain wall 
pinning depending on the lattice geometry and the den- 
sity of topological defects. In the 2d Potts model 
with g > 3 such point defects are the convergence of 
three distinct phases borders [H, The asymptotic 
density of defects that depends on q (negligible for g < 4 
but nonzero for g > 4) was also related with glassy be- 
havior [l^, In the presence of thermal fluctuations, 
the excess energy due to too many domain borders dis- 
appears after a transient 

In a sequence of papers j7l.ll9l.[20|. we explored the time 
evolution of domain and hull-enclosed area distributions 
in bidimensional coarsening systems with scalar order pa- 
rameter. To study the non-conserved order parameter 
case, on the one hand we used a continuum description 
based on the Ginzburg-Landau equation for the scalar 
field and the AUen-Cahn law dA/dt = —A. On the other 
hand, we studied the kinetic Ising model on a square 
lattice. For the q ~ 2 Ising case with non-conserved or- 
der parameter, the absence of triple points leads to the 
uniform shrinkage of all hull-enclosed areas. More con- 
cretely, the number of hull-enclosed areas per unit system 
area, in the interval {A, A + dA) at time t is related to 
the distribution at the initial time through 



nh{A,t) [A + \{t-ti),U 
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Once the initial distribution is known, the above equa- 
tion gives the distribution at any time t. For example, 
if one takes as initial states at = configurations in 
equilibrium state at Tc, where the distribution is exactly 
known for g = 2 [2l|, nh{A, 0) = Ch/A^ with ch = SttVS, 
the distribution at t > is nh[A,t) — Ch/{A + Xht)^, 
that compares extremely well with simulation data for 
the Ising model on a square lattice Q- For a quench 
from infinite temperature, on the other hand, the initial 
distribution corresponds to the one of the critical ran- 
dom continuous percolation 0, [l^ . This observation is 
an essential ingredient to understand the fact that, and 
obtain the probability with which, the system attains a 
striped frozen state at zero temperature (see 22| and 
references therein). Interestingly, on a square lattice, the 
initial state does not correspond to the random percola- 
tion critical point but the coarsening evolution gets very 
close to it after one or two Monte Carlo steps. There- 
fore, although Eq. ([3]) was obtained with the continuous 
description and, a priori, it is not guaranteed to apply 
on a lattice, it does describe the coarsening dynamics 
of the discrete Ising model remarkably accurately. In 
the presence of vertices, instead, the single hull-enclosed 
area evolution depends on n and, notably, a given hull- 
enclosed area can either shrink or grow depending on 
whether n < 6 or n > 6, respectively. Therefore, one 



cannot write a simple relation as the one in Eq. ([3]) to 
link the area distribution at time t to the one at the initial 
time ti and the distribution might get scrambled in a non- 
trivial way during the coarsening process (for example, 
when a domain disappears, the number of sides of the 
neighboring domains changes, along with their growth 
rate). 

In this paper we study the distribution of hull-enclosed 
areas during evolution in a Potts model with different 
number of states, notably g < 4 and g > 4. We then try 
to give an answer to some questions that can be posed 
when dealing with more than two competing ground 
states. To what extent the results obtained for g = 2 
are also valid for g > 3? In particular, in a quench from 
infinite temperature, does the random percolation critical 
point affect the evolution as it does for q — 27 Moreover, 
which is the interplay between nucleation and growth 
when the system goes through a first-order phase transi- 
tion as in the cases g > 4? What happens when the start- 
ing point of the quench is a first order transition point, 
with finite range correlations? In the presence of weak 
disorder, when first order transitions become continuous, 
are the scaling functions universal (super-universality)? 

The paper is organized as follows. In section [TT] we in- 
troduce the Potts model and, in section Hill we study the 
space-time correlation function starting from differently 
correlated initial states. Then, in section HVl several area 
distribution functions, either in the initial state or after 
the quench in temperature, are studied and described in 
detail. Afterwards, in section |Vl the Potts model with 
ferromagnetic random bonds (weak disorder) is studied 
and compared with the pure model. Finally, we make 
some final remarks and conclude. 



II. THE POTTS MODEL 



We consider the Potts model [15| on a square lattice of 
N = L'^ spins with L ~ 10'^ — 5 x 10'^, the state variables 
of which, Si with i = 1 ... TV, assume integer values from 
1 to g. The Hamiltonian is given by 
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(4) 



where the sum is over nearest-neighbor spins on the lat- 
tice and, until section |Vl Jij — 1, Vi,j. The transition, 
discontinuous for g > 4 and continuous for g < 4, oc- 
curs at Tc = 2/ln(l + y/q). We run from 500 to 4000 
Swendsen-Wang (SW) algorithm steps to reach an equi- 
librium initial condition at the critical point, and av- 
erage over 1000 samples to build the correlations and 
distributions. After equilibrium is attained and the sys- 
tem quenched, the evolution follows the heat-bath Monte 
Carlo algorithm [23| and all times are given in Monte 
Carlo step (MCs) units, each one corresponding to a 
sweep over N randomly chosen spins. 

After a quench in temperature, the abrupt increase 
in the local correlation creates an increasing order, de- 
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FIG. 1: (Color online) Rescaled space-time correlation func- 
tion at several times (from f = 2* to 2^^ MCs) after a quench 
from To — > oo to Ty = Tc/2 for q = 2, 3 and 8, with and with- 
out weak disorder (the disordered case is discussed in section 
rV)) . A good collapse is obtained for all q when rescaling the 
spatial variable by R{t) obtained from C(R,i) = 0.3 (note 
that the working temperature is different in all cases). Inset: 
R^{t) against t in a double logarithmic scale for q — 2, 3 and 
8 (from top to bottom). The characteristic length, related 
to the average domain radius, depends weakly on q and T 
for the pure model (through the prefactor). Within the time 
window explored there are still some small deviations from 
the expected AUen-Cahn behavior, R{t) ~ t^^^ in the case 
q = 8. When weak disorder is introduced (not shown), the 
growth rate is greatly reduced and R'' deviates from the linear 
behavior. 



spite the competition between different coexisting sta- 
ble phases. In the cases 9 < 4 the paramagnetic state 
becomes unstable at Tc and after a quench to any non- 
zero sub-critical temperature the system orders locally 
and progressively in patches of each of the equilibrium 
phases; in other words, it undergoes domain growth. In 
the cases q > A one can distinguish three working tem- 
perature regimes: at very low Tf the system gets easily 
pinned; at higher Tj but below the spinodal temperature 
Ts at which the paramagnetic solution becomes unstable 
the system undergoes domain growth; above the spinodal 
Ts there is competition between nucleation and growth 
and coarsening. However, the latter regime is very hard 
to access since Tg is very close to Tc (e.g. T, ~ 0.95Tc 
for g = 96 [24]). In this work we focus on the dynamics 
in the intermediate coarsening regime. 

The detailed evolution of the system during the coars- 
ening dynamics depends on the correlations already 
present in the initial state, whether absent (Tq 00), 
short- (Tc < To < 00 for all q and also To = Tc for 
g > 4) or long-range (Tq = Tc for g < 4). In the lat- 
ter case there is already one spanning cluster at t = 0, 
since the thermodynamic transition also corresponds to 
a percolation transition in 2d. On the other hand, for 
short-range initial correlations, such spanning domains 



are either formed very fast (e.g. in the case q = 2 with 
To —^ 00), or not formed at all (or at least not within the 
timescales considered here). 

III. EQUAL-TIME CORRELATIONS 

The degree of correlation between spins is measured by 
the equal time correlation function 

(5) 

where the average is over all pairs of sites a distance r 
apart. Away from the critical temperature, correlations 
are short-range, and after a quench from Tq > Tc to 
Tf < Tc these initial correlations become irrelevant after 
a finite time and the system looses memory of the initial 
state. In this sense, equilibrium states at all tempera- 
tures above Tc are equivalent. 

Figure [1] shows the correlation function C(r, t) as a 
function of the rescaled distance, r/R{t), after a quench 
from infinite temperature, where the initial correlation 
is null, to a working temperature Tf — Tc/2 for q = 2, 
3 and 8. In the inset, we show the length scale R{t) 
computed as the distance r at which the correlation has 
decayed to 0.3 of its initial value, that is, C{R,t) — 0.3. 
The linear behavior R"^ ^ t is clearly seen for g = 2. For 
q > 2 there are still some deviations from the t^^^ law 
at early times but these disappear at longer times (see 
also [Til [1^). Deviations at short times are indeed ex- 
pected, due to pinning at low temperatures (note that 
Tf = Tc/2 decreases with q). An excellent collapse of 
the spatial correlation is observed, as expected from the 
dynamic scaling hypothesis, when r is rescaled by the 
length R{t). Moreover, the universal curves seem to co- 
incide for these values of q, regardless of the different 
temperatures after t he q uench, in accordance with previ- 
ous evidence [H, [H, l26j for the fact that the space-time 
correlation scaling function (and its Fourier transform, 
the structure factor) are insensitive to the details of the 
underlying Hamiltonian. Whether this apparent q inde- 
pendence is only approximate or exact, and in the latter 
case whether it remains valid for very large values of g , 
are still open issues [13, HI] . As for the Ising model [1^ , 
the small r behavior is linear for all q, in agreement with 
Porod's law 

On the other hand, when the starting point is an equi- 
librium state at Tc with long-range correlations (g < 4) 
obtained by performing a sufficient number of Swendsen- 
Wang steps ^23j|, the system keeps memory of the long- 
range correlations present in the critical initial state. The 
equal time equilibrium correlation function decays at the 
critical temperature as a power-law, C(r, 0) ~ ^2-^-17^ 
where both the coefficient and the exponent 77 depend on 
g. For example, in two dimensions, 77 = 1/4 for g = 2 
and ?7 = 4/15forg = 3 [l^. Figure [5] shows the rescaled 
correlation function at several instants after the quench 
for g = 2. Some remarks are in order. First, there is a 
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FIG. 2; (Color online) Rescaled space-time correlation func- 
tion at several times (t — A, . . . , 2^° MCs) after a quench from 
To = Tc to Tf = Tc/2 for the Ising g = 2 case in log-log 
scale (cfr. Figs. |3] and |4] for 9 = 3 and g = 8, respectively). 
A very good collapse is obtained when rescaling the spatial 
variable by R{t), obtained from C{R,t) = 0.7. At t — 0, the 
correlation decays as r~'' (77 = 1/4), while for t > partial 
memory of the initial state is preserved, as evidenced by the 
power-law tail. The inset presents C(r, t) for the same times 
used in the main figure, with no rescaling of r. The correla- 
tion at f = is also shown, and deviations from the power-law 
behavior are already present, due to the strong fluctuations 
in magnetization at Tc- 



very good collapse for length scales up to r ^ R{t), where 
R{t) is such that C{R,t) = 0.7. Second, deviations from 
the initial power-law occur both at short and long length 
scales. Due to the ever growing structures, correlations 
decrease very slowly for small r. Indeed, C almost follows 
the plateau close to unity up to a certain, time increas- 
ing, distance. For long distances, on the other hand, 
these deviations occur because the initial states can be 
magnetized (at Tc, magnetization goes as L~^/'^) and, 
as r increases and the spins decorrelate, C(r, t) attains 
a plateau at rn^it). Since, eventually, the system equili- 
brates, the longer the time, the larger the magnetization 
and, consequently, the higher the plateau. Although not 
done here, it is also possible to postpone the approach 
to equilibrium by choosing initial states with very small 
magnetization, as done by Humayun and Bray [3l|, who 
also introduced a correction factor to account for bound- 
ary effects due to the system size being much smaller 
than the correlation length, L ^ ^, at Tc- Essentially 
the same behavior is obtained for q = 'd with 77 ~ 4/15 
(shown in Fig. [S]) and <? = 4 with 77 = 1/2 (not shown), 
respectively. We obtain the same scaling function for 
different final temperatures (not shown), indicating that 
super-scaling holds with relation to temperature. 

For g > 4, differently from the previous cases, there are 
no long-range correlations at T = Tc, only finite-range 
ones. The actual correlation length at Tc can be ob- 
tained analytically j3^ , but we did not attempt to mea- 
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FIG. 3: (Color online) Rescaled space-time correlation func- 
tion at several times (i = 4, . . . , 2^" MCs) after a quench from 
To = Tc to Tf = Tc/2 for g = 3 in log-log scale (cfr. Figs. [2] 
and 2] for q — 2 and g = 8, respectively). The slope r~'^ with 
f] — 4/15 is shown as a guide to the eye. Inset: R^{t) for q — 2 
(top) and 3 (bottom), obtained from C{R,t) = 0.7 together 
with the law R? c^t. 
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FIG. 4: (Color online) Rescaled space-time correlation func- 
tion at several times {t = 2^", . . . , 2^* MCs) after a quench 
from To = Tc to Tf — Tc/2 for g = 8 in linear-linear scale 
(cfr. Figs. [2] and [3] for q — 2 and g = 3, respectively). 



sure it, since the correlator that we use does not remove 
power-like prefactors [s^. For quenches below the limit 
of stability of the paramagnetic high-temperature state, 
see Fig. m for q = 8, the finite correlation length at t = 
is washed out once the scaling regime is attained. As 
might have been expected, the scaling function is indis- 
tinguishable, within our numerical accuracy, from the one 
for systems with a continuous phase transition (g = 2, 3) 
quenched from infinite temperature, cfr. Fig.[T] The col- 
lapse is not as good at short times since it takes longer for 
the system to approach the scaling regime in the q ^ 8 
case (see the inset). 
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After some discussion it is now well-established that 
weak quenched disorder changes the order of the phase 
transition from first-order to second-order when q > 4. 
The critical exponent 77 depends on q although very 
weakly. The scaling function of the space-time corre- 
lation function after a quench from the critical point also 
depends on q [H, . 



IV. AREA DISTRIBUTIONS 

More insight on the growing correlation observed dur- 
ing the coarsening dynamics after the quench can be 
gained from the study of the different area distributions 
n{A,t). We consider two measurements: geometric do- 
mains and hull-enclosed areas. Geometric domains are 
defined as the set of contiguous spins in the same state. 
In principle, such domains may enclose smaller ones that, 
in turn, may also enclose others and so on. One can also 
consider the external perimeter of the geometric domains 
(the hull) and the whole area enclosed by it. Let us dis- 
cuss the distribution of these areas after a quench from 
higher to low temperatures. 



A. The initial states 

Prior to the quench, the system is prepared in an initial 
state having either zero, finite or infinite range correla- 
tion, corresponding to Tq = 00 (Vg), Tq — Tc {q > 4) 
and Tq = Tc (9 < 4), respectively. Infinite temperature 
states are created by randomly assigning, to each spin, 
a value between 1 and g, while to obtain an equilibrium 
state at Tc, the system is further let evolve during a suf- 
ficient number of Swendsen-Wang Monte Carlo steps at 
Tc (typically 500). 
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FIG. 5: (Color online) Equilibrium domain (blue triangles) 
and hull-enclosed area (red squares) distributions at To ^ cxa 
for several values of g, corresponding to the random percola- 
tion problem at 1/g occupation density. The data are shown 
in linear-log scale. 
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When the temperature is infinite, neighboring spins 
are uncorrelated, and the configurations can be mapped 
onto those of the random percolation model with equal 
occupation probability p — l/q. On the square lattice 
considered here, this state is not critical for all values of 
q (criticality would require that one of the species density 
be 0.59). The highest concentration of a single species, 
0.5, is obtained for q = 2. This corresponds to the con- 
tinuous random percolation threshold but on a lattice, 
however, it is critical percolation on the triangular case 
only. Still, for the q = 2 Ising model on a square lat- 
tice, as shown in Q, the proximity with the percolation 
critical point strongly affects the system's evolution. 

Figure [5] shows the geometric domains and hull- 
enclosed area distributions n{A, 0) — qnip{A, l/q), where 




{A,p) is the corresponding size distribution for ran- 
dom percolation with occupation probability p = l/q. 
The factor q comes from the fact that the q species 
equally contribute to the distribution, while in the per- 
colation problem, a fraction 1 — 1/qoi the sites is empty. 
Unfortunately, n^p is not known analytically for general 
values of p away from pc- For q = 2 (not shown, see [l9l|), 
the system is close to the critical percolation point, and 
the distributions are power-law, n{A, 0) ~ A'^'^ with 
T = 2 for hull-enclosed areas and t > 2 for domain ar- 
eas, up to a large area cut-off where they cross over to 
exponential decays. For (7 = 8, deviations from linearity 
are already perceptible and, moreover, the hull-enclosed 
and domain area distributions can be resolved. For suf- 
ficiently low p (but probably valid for all p < Pr) the 
distribution tail is exponential [s^, Is^, [H, [s 
and the simulation data is compatible with 
nrp{A,p) ~ A~^ exp[—f{p)A], where the exponent 9 de- 
pends on the dimension only {6 = 1 for d = 2). Since 
all q states are equally present, the domains are typically 
smaller the larger the value of q and the occupation de- 
pendent coefficient f{p) thus increases for increasing q 
or decreasing p. The data for the different distributions 
shown in Fig. [5] can only be resolved for large values of 
A and small values of q. These differences are due to 
domains embedded into larger ones, already present in 
the initial state, that may remain at long times after the 
quench. For larger values of q the distributions get closer 
to an exponential in their full range of variation and it be- 
comes harder and harder to distinguish hull-enclosed and 
domain area probability distribution functions (pdfs) . 



2. To = Tc and 2 < g < 4 

When the transition is second order and the system 
is equilibrated at the critical temperature Tc, the distri- 
bution of both domains and hull-enclosed areas follow a 
power law. For example, for q = 2 in d = 2, the hull- 



enclosed area distribution is given by nh{A, 0) 



.(2) 
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FIG. 6: (Color online) Equilibrium power-law distribution of 
hull-enclosed areas at Tc for several values of q, from 2 to 
8. The line is cj^^/A^, from Eq. Q. Taking cjf' as a fitting 
parameter one finds very small values obeying c^**' < c^"^' < 
c^^^ (see text). For q > 4, the transition is discontinuous 
but nonetheless the distribution follows the power law until 
a crossover length, beyond which it decays exponentially, as 
shown in the inset for q = 8, 10 and, for comparison, q — 2 
(the horizontal data). 

where cjf^ = l/Sny^- 0.0229 f2ll. Generally, the hull- 
enclosed area distribution for g = 2, 3 and 4, is found to 
be 

^^.(^0)= ^''^']'^"' (6) 

as can be seen in Fig. [6l We choose to use the q—l pref- 
actor instead of q for consistency with our previous work. 
This can be done, however, with a small modification of 
the constant, {q — l)ch = qc'y^- Thus, each spin species 
contributes with the same share to the total distribution. 
Notice that, unless for q = 2, the value of c^^^ is not 
known exactly. A rough estimate of the constants can 
be obtained by taking them as fit parameters; we obtain 
cf^ ~ 0.0203 and c^f^ ~ 0.0192. The fit for g = 2 yields 
~ 0.0227, that compares well with the exact value. 
The constants obey the inequality cj^'*^ < c^^'' < c^^K 

For geometric domains, a relation similar to Eq. ([6]) is 
obtained for the area distribution at Tc and 2 < g < 4. 
The exponent is slightly larger than 2, t — 379/187 for 
q = 2 [H, m, and seems to be independent of q 
for q = 3, 4. The coefficients, however, are not exactly 
known (numerically, for q — 2, it is close to c]^^ [l9|). 



L, the crossover may or may not be observed due to the 
weakness of the transition. For example, for q = 5, 6 and 
8, estimates [H, [45| of ^ are, 2512, 159 and 24, respec- 
tively. For q = 5, f is indeed larger than the system size 
L that we are using and the system behaves as it were 
critical. Figure [H] also presents data for the hull-enclosed 
area distributions in models with q = 6, 8 and 10, where 
the deviations from the power-law are very clear, and 
occurring at smaller values of A for increasing q, as ex- 
pected. Equation © is thus vaHd up to this crossover 
length. The inset of Fig. [5] shows that for areas that are 
larger than a certain value A*, the distribution decays 
exponentially for q > 4. Thus, the general form of the 
hull-enclosed area distribution, apparently valid for all 
values of q at Tc, is 

r..(A0)=(i^e--4 (7) 

where aq = for q < 4. 

The area distribution of geometric domains is very sim- 
ilar to the hull-enclosed area one. None of the coefficients 

(2) 

, are exactly known; the numerical results for q — 2 

described in [l^ suggest that c''^^ is very close to c^^^^ and 
the study of the distributions for g = 3, 4 suggests that 
the similarity between c^''' and cj^''' holds for general q as 
well. 



B. The coarsening regime 

In Fig. [7] we show snapshots of the system at t = 2^° 
MCs after the quench for different values of q and Tq. 
When the initial condition is equilibrium at Tc for 2 < 
q < 4, the system presents larger domains than when 
the initial conditions are random, since the latter has an 
exponential instead of a power-law size distribution. On 
the other hand, for g > 4, when the correlation length at 
Tc is finite, there is no visible difference in the snapshots. 
In these figures, with the exception of some very small 
thermally induced fluctuations, there is no domain fully 
embedded inside another one, at variance with the q — 2 
behavior IQj. Nonetheless, in Fig. [SJ one such rare em- 
bedded domain occurs, being created after a coalescence 
process among the neighbors. Thus, for q > 2, due to 
this shortage of embedded domains (neglecting the small 
thermal fluctuations), there is little difference between 
geometric domains and hull-enclosed area measures, im- 
plying that the parameters of the distributions should 
have similar values. 



3. To =Tc andq> 4 

When the transition is discontinuous, the power law 
exists up to a crossover length, typically of the order of 
the correlation length ^, where the distribution deviates 
and falls off faster. Depending on the values of ^ and 



1. 2 < q < 4 and To ^ oo (finite correlation length) 

Despite the existence of many similarities between the 
Potts model with g = 2, 3 and 4, there are also some 
fundamental differences. Besides all having a continuous 
transition, the equal time correlation function seems to 
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FIG. 7: (Color online) Snapshots at t = 2^° MCs after a 
quench ft'om To ^ oo and To = Tc, left and right, respectively 
to Tf = Tc/2, for g = 3 [(a) and (b)] and g = 7 [(c) and (d)]. 
Different colors correspond to different species. Notice that, 
for q = 7, some domains have just coalesced (for example, the 
large yellow domain in panel (c), an effect that is absent for 
q = 2 and becomes rarer as q increases. Figure |8] gives a more 
detailed view of such occurrence. A few thermal fluctuations 
are also visible as small dots inside the clusters. Note that the 
structures are typically larger on the right column (To — Tc) 
than on the left one (To — > oo) for q = 3. This is not the case 
for q — 7 where there are no big structures in the initial state 
since ^{q — 7, Tc) is finite. 



share the same universal scaling function and the related 
correlation length grows with the same power of time, 
t^/^, see Sect.ini Nevertheless, differently from the q = 2 
case, that presents a percolating domain with probability 
almost one as early as i = 2 after the quench 0, [l^ , the 
q > 2, To — > cx) initial condition is sufhciently far from 
critical percolation that the system remains, at least in 
the time window of our simulations, distant from the per- 
colation threshold (in spite of the largest domain steadily, 
but slowly, increasing with time). 

As the system evolves after the quench, the distri- 
bution keeps memory of the initial state, that corre- 
sponds to random percolation with occupation probabil- 
ity p = l/q. And, by not getting close to a critical point, 
the distributions do not become critical and, as a con- 
sequence, do not develop a power law tail, as illustrated 
in Fig. [5] for q = 3 and Tq — > oo. There is, however, an 
envelope that is a direct consequence of dynamical 
scaling, present also for other values of q. The scaling hy- 
pothesis requires that the hull-enclosed area distribution 
satisfies nh{A,t) = t~'^nh{A/t). As a consequence, the 
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FIG. 8: (Color online) Snapshots at several times after a 
quench from To ^ oo to Tf — Tc/2, for q — 3. Different 
colors correspond to the three different species. The snap- 
shots show the coalescence of two white domains inside the 
100 X 100 zoomed region. Notice also the rare occurrence of a 
domain fully embedded into a single domain at t = 800 MCs. 



envelope straight line shown in the figure should have 
a —2 declivity. To see this, consider two of the curves 
shown (corresponding to times i* and a-nd define A* 
as the location of the point that is tangent to the enve- 
lope. If we use this value to rescale the distributions, we 
obtain R^{tl)nl = R'^{tl)nl, where n* is n{A* / {t*)) . 
That is, nl/n\ = [Al/Alf. Taking the logarithm of 
both sides, as in the graph, this gives the —2 declivity. 
The tangent point can thus be an alternative way to ob- 
tain the characteristic length in these systems. 

Very similar distributions (not shown) are obtained for 
g > 4 and Tq = Tc, another case with only finite correla- 
tion length in the initial state. 



2- 2 < g < 4 and To = Tc (infinite correlation length) 



Figure \W\ shows the collapsed distribution of hull- 
enclosed areas after a quench from Tq — Tc to Tf — Tc/2. 
We do not show the related figure for the geometric do- 
mains (without any spanning cluster) since they are al- 
most indistinguishable, another indication that the pa- 
rameters appearing in both distributions differ by very 
little. The overall behavior is similar to the q = 2 case, 
in which the collapse of curves for different times onto a 
single universal function demonstrates the existence of a 
single length scale that, moreover, follows the AUen-Cahn 
growth law, R{t) ^ i^/^. Assuming now that the number 
of sides in the von Neumann equation ([2]) can be replaced 
by a constant mean, n — *■ (n) , and using Eq. ([6]) and the 
results in Refs. 0, E3|, the hull-enclosed area distribu- 
tion for (7 < 4, within this mean-field-like approximation 
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has a different rate, either positive or negative, depend- 
ing on its number of sides, the average, effective A is 
constant, as is the case for the Ising model q = 2. Simi- 
larly, only those that have A„ff < present a power law 
distribution. Interestingly, the case with q — S with fer- 
romagnetic disorder, to be discussed in section |Vl seems 
to be marginal, Aeff — 0. Thus, there seems to be a net 
difference between cases that give a power law pdf the 
ones that do not. The detailed implications, that prob- 
ably involve the knowledge of perimeter and number of 
sides pdfs are beyond the scope of this paper and we 
postpone their study to a future work. 



FIG. 9; (Color online) Hull-enclosed area distribution at sev- 
eral times (given in the key) after a quench from equilibrium 
at To ^ oo to Tf = Tc/2 in the q = 2> case. Analogous dis- 
tributions are obtained for g > 4 and To ^ oo (not shown). 
The declivity of the envelope is —2 as a consequence of the 
scaling obeyed by the distribution (see text). 
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FIG. 11: (Color online) Average rate of area change as a func- 
tion of time for several cases studied in the text. Although the 
von Neumann law predicts that each domain has a different 
rate depending on its number of sides, the average is constant 
in time and non-positive for cases with a hull-enclosed and 
domain area pdfs with power-law tails. 



FIG. 10: (Color online) Collapsed hull-enclosed area distri- 
butions at several times after a quench from equilibrium at 
To = Tc to Tf = Tc/2, for g = 3. The line is Eq. ^ with 
X\f ^ ~ 1.4. The points at A/t <^ 1 that deviate from the 
scaling function are due to thermal fluctuations, depicted as 
a continuous black line in the inset. 



becomes 



nh{A, t) 



(9 -IK 



(9) 



(?)- 



2 ' 



(8) 



(3) 

that fits very accurately the data using Xf^ ~ 1.4. The 
deviations present at small values of A/t are due to ther- 
mal fluctuations that are visible in the inset (see fToJ for 
details). The above mean-field-like approximation can be 
directly tested by measuring the average change in area, 
dA/dt, that is, the number of spins included or excluded 
in those domains that survived during a given time in- 
terval. Figure [TT] shows the results for several cases. Al- 
though the von Neumann law predicts that each domain 
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FIG. 12: (Color online) Hull-enclosed area distribution after 
a quench from To — Tc to T/ = Tc/2 in the q — 8 model at 
different times given in the key. 
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3. q > 4: (finite correlation length) 

For g > 4, no initial state is critical and the initial 
distributions are not power law. After the quench, the 
system keeps some memory of this initial state and the 
scaling state has a more complex universal function. Fig- 
ure [T^] exemplifies this behavior for q = 8 and Tq — Tc- 
Again, since the quench is to Tc/2 the small A region, 
almost time independent, is due to the large number of 
small clusters formed by thermal fluctuations. After this 
initial decaying region, the distribution increases with A 
(the width of such region increases with time), to then de- 
crease exponentially. As for a quench from Tq — > oo and 
q = 3, an envelope power law forms when several 

distributions for different values of t are considered, and 
this is a direct consequence of dynamical scaling. The 
precise analytic form of the scaling function is not known 
and has been a matter of debate for several decades not 
only for the Potts model but for related models of inter- 
est for the grain growth community (see [H, 113, IH, |49| 
and references therein), this issue being still unsettled. 

The tail keeps memory of the initial condition and 
is exponential for q > A. This fact is better appreci- 
ated in Fig. [T3] where we use a scaling description of 
the data in Fig. [12] and a linear-log scale. The univer- 
sal scaling function shown in the inset, for A > R^{t), 
is the same for both Tq — Tc and To — > cxd (not shown): 
f{x) ^ exp(— aa:), with a fitting parameter a ~ 0.23 
in both cases. This is a consequence of the finiteness of 
the correlation length at i = and the fact that the tail 
of the distribution samples large areas. Notice also that 
the collapse is worse for small values of A/i?^(i), prob- 
ably due to the temperature roughening of the domain 
walls. 



V. WEAK QUENCHED DISORDER 

It is also possible to study the coarsening behavior in 
the presence of random ferromagnetic (weak) disorder ob- 
tained, for example, by choosing the bonds from a prob- 
ability distribution function, P(Jy ), with semi-definite 
positive support, > 0. Weak disorder weakens the 
phase transition that, in the 2d random bond Potts model 
(RBPM) becomes continuous for all g [13, [HI, [13, [H [Hi . 
For the bimodal distribution 



1 



1 



P{J^j) = - Ji) + T^K'h, - J2), 



that we shall treat here, the transition occurs at [5£ 



(9) 



1) (e^-'- -l)=q. 



(10) 



We use Ji = 1 and J2 = 1/2, that is, J1/J2 = 2, for 
which Tc ~ 1.443 for g = 3 and T^ ~ 1.087 for g = 8. 
The additional contribution to interface pinning due to 
quenched disorder may be useful to understand phenom- 
ena such as the so-called Zener pinning [5^ . [stI [58| . 




A/R\t) 



FIG. 13: (Color online) Dynamical scaling of the distributions 
shown in Fig. [12] for g = 8 and To = Tc. A good collapse is 
obtained, for A > R^{t), with nh{A,t) = R-'^{t)f{A/R^) 
(the data for A < 10, corresponding to thermal fluctuations, 
have been removed from the inset). The main panel shows 
the exponential tail, f{x) ~ exp{—ax), where the fitting 
parameter is a ~ 0.23. Data for To <x collapse onto the 
same universal function (not shown) . The inset also shows the 
small area region, where the collapse is not as good, probably 
due to thermal interface roughening. 



The growth law for the random q = 2 case, whether 
logarithmic or a power-law with a T and disorder- 
strength dependent exponent, has been a subject of de- 
bate |59l . 160 . [giI . [6^ and arguments for a crossover be- 
tween the pure growth law t^^^ to logarithmic growth 
at an equilibrium length-scale - that can be easily con- 
fused with the existence of T and disorder dependent 
exponent in a power-law growth - were recently given 
in [g^ . Whether the asymptotic logarithmic growth also 
applies in the case q > 2 has not been tested numeri- 
cally. Regardless of the growth law, dynamical scaling 
is observed in the correlation and area distribution func- 
tions for q — 2 [g^I as well as for g > 2. The rescaled 
space-time correlations for the disordered q — 3 and q — 8 
after a quench from Tq ^ 00 are indistinguishable from 
the pure case, see Fig. [TJ Notice that the initial con- 
dition, being random, is not affected by the presence of 
disorder, differently from the To = Tc case, where the 
disorder changes the nature (correlations) of the initial 
states, leading to a failure of the super-universality hy- 
pothesis. In this case, long range correlations are present 
and the decay is power-law, C(r) ^ r~'', as shown in 
Fig. [13] in analogy with the pure 2 < q < 4 cases. How- 



65l [66j , for the bi- 



ever, as it was pointed out in Refs. 
modal distribution and the values of Ji chosen here, this 
may suffer from crossover effects, not yet being in the 
asymptotic scaling regime. For q — 3 and 8, the fitted 
exponent is 77 ~ 0.27 and 0.28, respectively, results that 
are compatible with Refs. (gsI. [67j. 

Figure [15] shows the rescaled hull-enclosed area distri- 
bution function for the disordered q = 8 model quenched 
from equilibrium at Tq = Tc- The transition is continu- 
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FIG. 14: (Color online) Rescaled space-time correlation func- 



tion at several times {t — 2 



2^° MCs) after a quench 
from To = Tc to Tf = Tc/2 for q = 8 and bimodal disorder 
in log-log scale (compare with the pure q = 8 case in Fig. |4]). 
The tail keeps memory of the r~'^ decay of the initial equi- 
librium state. We found that 77 ~ 0.28 fits well the data. A 
similar figure is obtained for q = 3, with an exponent close to 
0.27. Inset; equilibrium distribution of hull-enclosed areas at 
Tc with weak disorder (red) - a critical case with long-range 
correlations - and without (blue) - a case with short-range 
correlations. The initial distribution was obtained after 3000 
SW steps and is indistinguishable, within the errors, from the 
distribution of the pure model with g = 3 (line) . 



ous for all q and the equilibrium area distribution at Tc 
is described, for large A and within the simulation error, 
by the same power law as the pure model with q — 3, 
Eq. dS]), see the inset of Fig. [TH There are, however, 
small deviations for not so large values of A that, after 
the quench, become more apparent, as can be observed in 
the main panel of Fig. [121 ^^^d in the inset where a zoom 
on this region is shown. This may be an effect of pinning 
by disorder that keeps an excess of small domains, slow- 
ing down their evolution, and not letting dynamic scaling 
establish at these scales. The weak increase of the pdf at 
small As, where the type of dynamics is important, re- 
sembles also the behavior of the pure case, see the inset 
in Fig. [131 

Interestingly enough, the long-range correlations in the 
Tc initial condition, that in presence of quenched disorder 
exist for all q, determine the large area behavior of the 
dynamic pdf, more precisely, they dictate their power-law 
decay. Thus, for q > 4 pure and disordered systems be- 
have qualitatively differently when the initial conditions 
are in equilibrium at Tq = Tc. An interesting character- 
istic of both figs.linifor q — % and [TBI for q — ?> \s that the 
tail is well described by the pure distribution with g = 3. 

For g < 4, the initial state at Tc has long-range corre- 
lations with and without disorder. We study the effect of 
disorder on the hull-enclosed area distribution in Fig. 1161 
The tail is a power law compatible with the —2 expec- 
tation. The small area behavior is interesting: the data 



deviates from the critical behavior as can be seen in the 
inset of Fig. 1151 showing a non-monotonic excess of small 
domains that resembles the behavior of the pure case. 
Somehow, disorder, even if local, has a strong effect at 
large scales while affecting less the properties at small 
ones. 




FIG. 15: (Color online) Collapsed hull-enclosed area distribu- 
tion for q = 8 after a quench from To = Tc to Tj — Tc/2 in 
the RBPM (half the bonds are changed from 1 to 0.5). Data 
taken at several times are shown with different symbols. The 
thick line shown is 20^^^ /{I + x)^ (the pure q = 2> case). The 
tail keeps memory of the initial state and is well described 
by a power law (solid red line) with exponent —2. Inset: the 
small A region does not scale and its behavior is reminiscent 
of what is seen in the pure case at small areas (but narrower), 
see the inset in Fig. 1131 



VI. CONCLUSIONS 

We presented a systematic study of some geometric 
properties of the Potts model during the coarsening dy- 
namics after a sudden quench in temperature. Although 
the theory for the Ising model (g = 2) cannot be easily 
extended to g > 2, our numerical results are a first step 
in this direction. 

Our analysis demonstrates the fundamental role played 
by the initial conditions, more precisely, whether they 
have an infinite correlation length or not. 

The distribution of hull-enclosed areas in pure 2d Potts 
models with 2 < g < 4 and dirty cases with all g that 
evolved from initial conditions in equilibrium at Tc, that 
is to say cases with an infinite correlation length, are 
large, with power-law tails that are compatible with the 
exponent —2 within our numerical accuracy. The small 
area behavior is richer. In cases with 2 < g < 4, again 
within our numerical accuracy, the hull-enclosed area is 
well captured by a simple extension of the distribution 
in Eq. ([U]): 



nh{A,t) 



(9 



[A + R^{t)Y 



2<'Z<4, (11) 
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FIG. 16: (Color online) Collapsed hull-enclosed area distri- 
bution for q = 3 after a quench from To = Tc to T/ = Tc/2 
in the RBPM (half the bonds are changed from 1 to 0.5). 
Differently from the q = 8 case shown in Fig. 1151 here both 
the pure and disordered initial state are critical (inset), and 
the distributions at f = are indistinguishable. Nonetheless, 
the subsequent evolution of these cases are not equivalent. 
For roughly A < Ft?{t) the scaling function for the disordered 
case no longer follows the pure case function, 20^^^ /{I + x)'^ . 
Interestingly, this equation describes the tail for both q = 3 
and 8 (previous figure). 

that upgrades the prefactor Ch to depend on q and dis- 
order. In disordered cases with g > 4 this form does not 
describe the small A dependence. Indeed, the scaled pdf 
has a non monotonic behavior both with and without dis- 
order, with dynamic scaUng faihng in the explored times. 



possibly due to strong pinning effects. 

Very different are the distributions of dynamic hull- 
enclosed areas evolved from initial states with finite cor- 
relation lengths, as those obtained in equilibrium at 
To — *■ cxD. The scaling functions of the dynamic distri- 
butions are reminiscent of the disordered state at the 
initial temperature with an exponential tail. Differently 
from the Ising case, in which the percolation critical point 
gave a long-tail to the low-temperature distribution, in 
cases with q> 2 this does not occur. 

There are further geometrical properties that were not 
explored here and deserve attention. For example, the 
distribution of perimeter lengths, number of sides, and 
the correlation between area and number of sides are of 
interest in cellular systems. In particular, it would be 
interesting to check whether the Aboav-Weaire and Lewis 
law [68] are valid in the 2d Potts model with and without 
random bonds. 

Finally, the Potts model is realized in a few experi- 
mental situations. Our results should be a guideline to 
search for similar distributions, in analogy to what has 
been done in the q = 2 case, by using a liquid crystal 
sample |i69,] . 
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